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ABSTRACT 

We present sensitive high angular resolution (0.57" -0.78") SO, S0 2 , CO, C 2 H 5 OH, HC 3 N, 
and HCOCH2OH line observations at millimeter and submillimeter wavelengths of the young 
O-type protostar W51 North made with the Submillimeter Array (SMA). We report the presence 
of a large (of about 8000 AU) and hot molecular circumstellar disk around this object, which 
connects the inner dusty disk with the molecular ring or toroid reported recently, and confirms 
the existence of a single bipolar outflow emanating from this object. The molecular emission 
from the large disk is observed in layers with the transitions characterized by high excitation 
temperatures in their lower energy states (up to 1512 K) being concentrated closer to the central 
massive protostar. The molecular emission from those transitions with low or moderate excita- 
tion temperatures are found in the outermost parts of the disk and exhibits an inner cavity with 
an angular size of around 0.7". We modeled all lines with a Local Thermodynamic Equilibrium 
(LTE) synthetic spectra. A detail study of the kinematics of the molecular gas together with a 
LTE model of a circumstellar disk shows that the innermost parts of the disk are also Keplerian 
plus a contracting velocity. The emission of the HCOCH2OH reveals the possible presence of 
a warm "companion" located to the northeast of the disk, however its nature is unclear. The 
emission of the SO and SO2 is observed in the circumstellar disk as well as in the outflow. We 
suggest that the massive protostar W5 1 North appears to be in a phase before the presence of 
a Hypercompact or an Ultracompact HII (HC/UCHII) region, and propose a possible sequence 
on the formation of the massive stars. 
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1. Introduction 



It is thought that the Hypercompact or Ultracompact HII (HCHII or UCHII) regions trace early stages 
of massive star formation after the hot molecular core phas e. However , the physical p rocesses involved in 
this transition are still poorly understood (see for reviews: lKurtzll2005l : iLizanol 120081) . The appearance of 
one HCHII begins when the Lyman continuum output from a massive young star becomes sufficient high 
to ionize its surroundings. These objects show a rising continuum spectra, S v ocv a , with a slope a ~ 1, 
intermediate between the optically thick and thin limits, small sizes (~ 0.01 pc), high temperatures (10 4 K), 
high densities (> 10 cm ), and very broad radio recombination lines (Avfw hm > 40 km s" 1 ) i ndicating 
both pressure broadening and the presence of bulk motions of dense gas (e.g. ISewilo et al.ll20041) . On the 
other hand, the hot molecular cores internally heated by massive protostars are small (~ 0.05 pc) and compact 
regi ons of dense molecular gas (> 10 6 cm" 3 ) which sh ow very rich millimeter and submillimeter spectra 
(see lKurtz et al.ll2000l : Ivan der Takll2004l : ICesaronill2005L for reviews on this topic). The hot molecular core 
phase begins at much colder temperatures (~ 100 K) than the HC/UCHII stage. 

There are many cases where the UCHI I or even the H CHII regio ns are observed in a ssociation with 
hot cores, so me of them are: G5 . 89-0.39 Jsu et al.l l2009h: G24 Al teeltran et al.1 l2006h : G10.47, and 
G31.41+0.31 Jcesaroni et aDboidlRo lffs et alJboOSlh: W 51e2. W51e8, NGC 7538 IRS1, G28.2 and G10.6 



( Klaassen et al 



aoaJTlRAS 17233-3606 dLeurini et al.ll2008h . However, there are some others cases where 



the hot cores do not show the presence of a strong HCHII (e.g. AFGL490, NGC 7538S, and IRAS 20126; 



Nakamura et al.lll99ll : ISchreyer et al.ll2002l : ICesaroni et al.ll2005l :ISandell & Wright 



Walmslev 



1995; 



2010 ) possibly beca use 



Osorioetal.il 19991 ) or 



the infalling molecular gas partially quenches its formation (e.g. 
maybe because the massive central protostar(s) are/is in an early phase with not high enough temperatures 
to ionize their/its surroundings. 

W51 North is one of the youngest massive stars within the luminous cluster W51-IRS2 that is still on 
the process of formation. This object is asso ciated with a stro ng mm and submm source, a hot molecular 
core at an excitation temperature of ~ 200 K (|Zhang et al.lll998r) . no centi meter free-free emission at all (< 



5 mJy at 1.3 cm) that co uld be related with an HC/UCHII or a thermal jet (|Gaume et al.ll 19931: IZhang et al. 



1998: Eisner et al.ll2002r) . From this object emanates a powerful molecular outflow obse rved at very small 



scales by masers s pots of SiO and FLO, and at large s cales in th ermal SiO(7 = 5-4 ) dSchneps et al.lll981 



Eisner et all2002l ; llmai et~al]|2002l ; Eapata et all2009l) . Recently. Eapata et al.l d2009h using high angular 1.3 
and 0.7 mm continuum, and SO2 line observations made with the Submillimeter Array (SMA) and the Very 
Large Array (VLA) resolved the large d usty and molecular hot core in a rotating Keplerian ring with a size 
of ~ 8000 AU (at a distance of ~ 6 kpc: llmai et ZM2002l : liarbosa et alEoollxu et al.ll2009h . This structure 
surrounds a compact central dusty circumstellar disk with a smaller size of approximately 3000 AU, and 
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where the molecular northwest(redshifted)-southeast(blueshifted) bipolar outflow emerges. However, the 
relationship between both circumstellar structures was not clear. 

In this paper, we present sensitive high angular resolution SMA submillimeter and millimeter line 
observations of the W5 1 North region that were made in an attempt to understand the nature of the circum- 
stellar molecular and dusty structures associated with this young massive protostar. In Section 2 discuss 
the observations undertook in this study. In Section 3, we present and discuss our SMA millimeter and 
submillimeter data. Finally, in Section 4, we give the main conclusions of the observations presented here. 

Table 1 : Observational and physical parameters of the submillimeter and millimeter lines 





Rest frequency 


^ lower 


Range of Velocities 


Linewidth 


LSR Velocity^ 


Peak Flux 


Lines 


[GHz] 


[K] 


Perns" 1 ] 


[kms- 1 ] 


[km s~ ! ] 


[Jy Beam ] 


C 2 H 5 OH(37 8 ,29-36<), 2 8) 


216.8121... 


715 


+54,+67 


7 


59 


0.1 


HCOCH 2 OH(62 13 ,4 9 -62 12 ,5 ) 


216.7585... 


1476 


+54,+66 


8 


59 


0.4 


CO(J = 3 - 2£ 


345.7959... 


16 


+20, +120 


50 


60 


3.0 


SO(8 9 -7 8 £ 


346.5284... 


21 


+30.+90 


35 


59 


1.0 


S0 2 (19i,i9-18o.is) 


346.6521... 


151 


+41.+72 


15 


58 


1.0 


HC 3 N(38-37) (v,=0) 


345.6090... 


306 


+51, +72 


9 


60 


2.0 


HCOCH 2 OH(68 19,49 -67 20 ,4 8 ) 


337.3988... 


1512 


+48,+73 


10 


62 


2.0 



"The linewidth and LSR velocity were obtained fitting a Gaussian to the spectra. 
*Spectra with multiple peaks and in absorption about a LSR velocity of 60 km s _I 
^Spectra with two peaks. 



2. Observations 
2.1. Millimeter 

The observations were obtained with the SMA0 during 2008 January 17. The SMA was in its very 
extended configuration, which included 28 independent baselines ranging in projected length from 30 to 
385 k/l. The phase reference center of the observations was R.A. = 19h23m40.05s, decl.= 14°31'05.0" 
(J2000.0). The frequency was centered at 217.1049 GHz in the Lower Sideband (LSB), while the Upper 
Sideband (USB) was centered at 228.1049 GHz. The primary beam of the SMA at 230 GHz has a FWHM 
of about 60". 

We detected the lines C 2 H50H(37 8j 29-36 9 ,28) and HCOCH20H(62i 3i 49-62i2,5o) in the LSB. See Table 
1 for their rest frequencies and rotational temperatures above of the ground state. The full bandwidth of 
the SMA correlator was 4 GHz (2 GHz in each band). The SMA digital correlator was configured in 24 
spectral windows ("chunks") of 104 MHz each, with 256 channels distributed over each spectral window, 



'The Submillimeter Array (SMA) is a joint project between the Smithsonian Astrophysical Observatory and the Academia 
Sinica Institute of Astronomy and Astrophysics, and is funded by the Smithsonian Institution and the Academia Sinica. 



-4- 



providing a resolution of 0.40 MHz (0.56 km s *) per channel. However, in this study we smoothed our 
spectral resolution to about 1 km s" 1 . 

The zenith opacity (t 2 3og#z)> measured with the NRAO tipping radiometer located at the Caltech Sub- 
millimeter Observatory, was ~ 0.15, indicating good weather conditions during the observations. Observa- 
tions of Uranus provided the absolute scale for the flux density calibration. Phase and amplitude calibrators 
were the quasars 1925+211 and 2035+109. The uncertainty in the flux scale is estimated to be 15-20%, 
based on the SMA monitoring of quasars. Observations of Uranus provided the absolute scale for the flux 
density calibra tion. Further technical descriptions of the SMA and its calibration schemes can found in 



Ho et all (2004). 



The data were calibrated us ing the IDL superset MIR, originally developed for the Owens Valley Radio 
Observatory (IScoville et al Il993h and adapted for the SMA@The calibrated data were imaged and analyzed 



in the standard manner using the MIRIAD and AIPS packages. We used the ROBUST parameter of the 
INVERT task of MIRIAD set to 2 to obtain a slightly better sensitivity while losing some angular resolution. 
The resulting image rms noise of line images was around 30 mJy beam" 1 for each velocity channel (with a 
smoothed size of 1 km s" 1 ) at an angular resolution of 0"57 X 0"42 with a P.A. = 57.6°. 

Table 2: Observational parameters of the circumstellar multilayer disk 





Central Position 


Flux Density 


Deconvolved Size" 


Approx. cavity size 




or(2000) 


5(2000) 










Specie 


[h m s] 


[0 " '] 


[Jy Beam - ' km s ] 


[arcsec 2 ] 


[Degrees] 


[arcsec 2 ] 


SO(8,-7 8 ) 


19 23 40.074 


14 31 05.70 


341 ±20 


1.78±0.01 x 1.39*0.02 


171+2 


-0.8 


SO 2 (19i,i9-18 ,i8) 


19 23 40.071 


14 31 05.54 


234±15 


1.65±0.04 xl.49±0.03 


24±6 


-0.7 


HC 3 N(38-37) (v,=0) 


19 23 40.060 


14 31 05.54 


88±15 


1.36±0.04 x0.93±0.04 


135+7 


- 0.6 


C 2 H 5 OH(37 8 , 2 9-3 6 9 ,28) 


19 23 40.041 


14 31 05.49 


116+15 


1.31+0.03 x0.89±0.03 


107±7 


-0.8 


HCOCH 2 OH(68 19,49 -67 20 ,48 ) 


19 23 40.054 


14 31 05.44 


105+13 


1.14±0.03 x0.88±0.03 


99±3 




HCOCH 2 OH(62i3,49-62i 2j5 o) 


19 23 40.062 


14 31 05.59 


9.1±2 


0.81±0.03 x 0.40±0.09 


75±3 




— Companion 


19 23 40.077 


14 31 05.63 


8.5±1 


1.27±0.04 x0.39±0.03 


74±2 










Dusty circumstellar disk^ 








7 mm 


19 23 40.057 


14 31 05.67 




0.58±0.02 x0.27±0.02 


7()±6 





The Gaussian fitting was obtained with the task JMFIT of AIPS. 
'Data obtained from Zap ata et al. (2009J. 



2.2. Submillimeter 



The observations were obtained with the SMA on 2008 July 13. At the time of these observations the 
SMA had seven antennas in its extended configuration with baselines ranging in projected length from 30 



2 The MIR-IDL cookbook by C. Qi can be found at http://cfa-www.harvard.edu/~cqi/mircook.html 



-5- 



to 255 kA. The primary beam of the SMA at 340 GHz has a FWHM of 37". The molecular emission from 
the hot core was found well inside of the primary beam. 

The receivers were tuned to a frequency of 345.796 GHz in the upper sideband (USB), while the lower 
sideband (LSB) was centered on 335.796 GHz. The LSB contained line HCOCH 2 OH(68i 9j 49-672o,48) while 
the USB the lines SO(8 9 -7 8 ), S0 2 (19 U9 -18 a i8), CO(/ - 3 - 2), and HC 3 N(38-37) (v ; =0). See Table 1 for 
their rest frequencies and lower level energies above the ground energy state. The SMA digital correlator 
was configured in 24 spectral windows ("chunks") of 104 MHz each, with 128 channels distributed over each 
spectral window, thus providing a spectral resolution of 0.81 MHz (0.70 km s" 1 ) per channel. However, we 
smoothed our spectral resolution to 1.0 km s _1 per spectral channel. 

The zenith opacity (T230GHz) was ~ 0.035 - 0.04, indicating excellent weather conditions. Observations 
of titan provided the absolute scale for the flux density calibration. Phase and amplitude calibrators were the 
quasars 1751+096, and 1925+211. The uncertainty in the flux scale is also estimated to be 15 - 20%, also 
based on the SMA monitoring of the quasars. 

The calibrated data were imaged and analyzed in standard manner using the MIRIAD, and AIPS pack- 
ages. We also used the ROBUST parameter set to 2. The line image rms noise was around 100 mJy beam -1 
for each velocity channel (with a smoothed size of 1 km s~ J ) at an angular resolution of 0"78 x 0"58 with 
a P. A. - 84.7°. 



3. Results 



In Figure 1, we present the disk/ring/outflow system found tow ards W5 1 No r th usi ng observations of 
the SiO(/ = 5-4), 802(222,20-221,21), and millimeter continuum by IZapata et al. 1 i2009h . However, in this 



image, we have exchanged the SiO emission by the emission of the CO(/ = 3-2) presented here for the first 
time. The CO(7 = 3-2) confirms the presence of a single powerful and collimated (15°) bipolar outflow 
(with the blueshifted side towards the southeast while the redshifted one to the north) emanating from the 
disk/ring structure at similar scales. However, the CO(/ = 3-2) reveals that the outflow extents for a few 
arcseconds more on both sides in comparison with the emission traced by SiO(/ = 5-4), and shows a 
slightly larger radial velocity range (+30 to +120 km s" 1 ). In this image, it is more clear to see the deviation 



of the blueshifted side of the outflow with respect to the redshifted one than in lZapata et all (|2009). 



Figure 2 shows the spectra of the multiple lines reported here with exception of th e CO(7 = 3 - 2) . 
Additionally, we present synthetic LTE models for the spectra using XCLASS program (|Comito et al 1I2005I) . 



The data were fitted well with the synthetic model using rotational temperatures, column densities, and sizes 
of 250 K, 2 x 10 16 cm 2 , and 1", respectively, for the lines SO(8 9 -7 8 ), SO 2 (19 U9 -18 ,i 8 ), HC 3 N(38-37) 
(v f =0), and C 2 H 5 OH(37 8i29 -36 9 ,28) and 800 K, 5 x 10 18 cm 2 , and 0.5" for transitions with high excitation 
temperatures in the lower energy states, HCOCH20H(68i 9 ,4 9 -672o,4 8 ) and HCOCH20H(62i34 9 -62i2,5o)- 
This modeling serves as for avoiding the line confusion, and suggest that hot gas is only found in the 
innermost parts of the hot core. See Table 1 for the rest frequencies and the physical parameters of the lines. 
The physical parameters of all lines were obtained fitting a Gaussian to the line profile. 
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In order to be more confident about the detection of the glycol-aldehyde (HCOCH2OH) toward W5 1 
North, we have fitted simultaneously many other lines that fall in our 2 GHz bands with our synthetic LTE 
model using the same parameters as above. We obtained a reasonable agreement between the synthetic 
model and the detected lines. A full astrochemical analysis of these lines (and other species) is beyond the 
scope of the present paper in which we concentrate on the spatial distributions of the lines presented in Table 
1. 

The spectra of the SO(8 g-7s) additionally s hows self-absorption possibly due to the molecular gas is 
falling into the central object (IZapata et al. 2008). 



The integrated intensity or moment zero maps of the lines associated with the molecular disk are pre- 
sented in Figure 3. We have divided these lines into three groups showing the transitions with energy levels 
with moderate (20-160 K), high (300-800 K), and very high (1400-1600 K) energies and presented them 
in two different panels. From these images, it is easy to see how the the molecular emission is found in 
layers with the transitions characterized by high excitation temperatures in their lower energy states (up to 
1512 K), and/or critical densities, being concentrated closer to the central high-mass protostar W51 North. 
This phenomenon confirms clearly that the heating in W5 1 North is internal. Additionally, molecules with 
transitions characterized by low and moderate excitation temperatures i.e. SO(89-7g), S02(19i,i9-18o,is), 
HC3N(38-37) (v f =0), and C2H50H(37g,29-369,28) show a central cavity of approximately an angular size of 
0.7", see Table 2. These lines t race thus the molecular Keplerian ring (or toroid) around W51 North reported 
in SO2 by IZapata et all (|2009i) . see Figure 1. The transition of the molecule SO2 utilized by these authors 
has a similar excitation temperature in the lower energy state as those mentioned above. Note the similarity 
between the angula r size of the inner cavity and the size of the compact dusty circumstellar disk reported by 
Zapata et all d2009h . see Table 2. 



The molecular emission from the transition with a high excitation temperature in the lower energy 
states (i.e. HCOCH20H(68 19,49-6720,48)) is found peaking at the position of the cavity and showing similar 
angular sizes to this, see Table 2. This molecule is thus more intimately related with the circumstellar 
compact dusty disk. 

The observational parameters of all lines are shown in Table 2. We noted from Table 1 and 2, there 
is a correlation between the deconvolved sizes of the molecular emission and the excitation temperatures 
in lower energy states of each molecular specie. The molecular emission from transitions characterized by 
high excitation temperatures show to be very compact. This correlation could be obviously obtained for 
lines that are optically thick with T& ~ rj T ex , where the T^ is the brightness temperature, 77 is the filling 
factor, and T ex is the excitation temperature. 

In the left panel of Figure 3 the emission of the lines HCOCH 2 OH(62i3,49-62i 2 ,5o) and HC 3 N(38-37) 
(Vf=0) reveals the presence of a possible warm "companion" located to the northeast of the disk. We have 
marked the position of this putative com panion with a yellow triangle. This "companion" is also observed in 



SO2 a t the same position, see Figure 3 of IZapata et all (120091) . and is associated with a group of water maser 
spots ( Imai et al.ll2002 : Eisner et aDl2002h . The position of the "companion" was found by fitting a Guassian 



to its HCOCH20H(62i3,49-62i2,5o)) emission. However, it still not clear if the "companion" is real protostar 
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or if this could be the result of the interaction of the outflow with a high density zone of the molecular cloud. 



We do not find any clea r evidence of outflowing gas activity associated with this possible "companion" 
(see also lZapata et al.ll2009r) . Although the blueshifted side of the CO(/ = 3-2) outflow has not the same 
position angle as the redshifted one, this seems not be ejected from the companion. We show how the 
blueshifted side appears to be more likely ejected from W51 North in Figure 1. We drew a line that crosses 
this side of the outflow and points directly to the dusty compact disk. The SO2 and SiO emission indeed 
show how the outflow is deviated to where the CO(/ = 3 - 2) is located ( Zapata et all2009 ). However, more 
observations in some other molecular outflow tracers are thus necessity to firmly discard the existence of a 
second outflow in W5 1 North. It is interesting to note that the possible companion is not observed in the 
hotter molecular gas tracer HCOCH20H(68 19,49-6720,48) suggesting that it may not as warm as the central 
massive star. 

Figures 4 and 5 show the position-velocity diagrams (PV-diagrams) of different molecules which trace 
distinct scales of the disk. In Figure 3, we have marked the orientation and position of the PV-cuts. In the 
left panel of Figure 3 we show a white line with a PA. =30° that corresponds to the PV-cuts shown in Figure 
4. On the other hand, the white line with a PA.=160° in the right panel corresponds to the the PV-cuts 
shown in Figure 5. 

Figure 4 shows the molecular emission from HC3N(38-37) and HCOCH20H(68 19,49-6720,48) located 
in the innermost parts of the disk, while the PV-diagrams of molecules as SO(89-7g), and S02(19ij9-18o,i8)> 
which trace the outermost parts, are presented in Figure 5. The PV-diagrams in Figure 4 reveal that the hot 
molecular gas closer to the massive protostar i s Keplerian. In addit ion, in this figure we have overlaid the PV- 
diagram of the LTE Keplerian disk modeled in lZapata et al.l (120091) . but without an inner cavity and a smaller 
size. Both structures shown a very good correspondence. The SO(8 9-7s) and S02(19i.i 9-18q.is) presented 
in Figure 5 trace much larger structures similar to the ring reported in lZapata et al.1 (|2009). These molecules 
in addition show clearly two northwest and southeast high velocity extensions excited by the bipolar outflow 
mapped in CO(7 = 3-2) and SiO(7 = 5-4). The PV-diagram of the molecule C 2 H 5 OH(37 8j2 9-369 ; 28) 
shows also Keplerian motions, while that obtained from the molecule HCOCH20H(62i3 49-62i2,5o) shows 
a more compact structure without a clear morphology, these diagrams are not presented in this study. 

Finally, in Figure 6 we show the spectral energy distribution (SE P) from the c entim eter, millimeter, 
and submillimeter wavelengths of W51 North with data obtained from IZapata et al.l (120091) and Tang et al. 
(in prep.). The new value at submillimeter wavelengths presen ted here shows also to be good fitted with the 
a = 2.8 value (where the flux density goes as S v -v a ) found by lZapata et all (|2009h . 
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4. Discussion 

4.1. A multilayer infalling and keplerian disk around W51 North: Connecting the molecular ring 

and dusty compact disk. 



The combined millimeter and s ubmillimeter l ine ob servations presented here from the ring and com- 
pact dusty disk structures reported in lZapata et all (120091) suggest that we are seeing a single hot multilayer 
accreting disk around W51 North (Figure 3). This large molecular disk connects both circumstellar struc- 
tures, the dusty compact disk, observed at 7 and 1 mm, and the molecular ring mapped in SO2 in a single 
structure. Some molecules are found in association with the dusty compact disk while some others with the 
ring. 

The molecular emission associated with the multilayer disk reveals that the molecular emission with 
transitions with a low excitation temperature in the lower energy states delineates the outer disk and shows 
a cavity, while such transitions with a high excitation temperature in the lower energy state are found only 
in the innermost parts of the disk and closer to the central massive protostar(s). One would expect to this 
physical phenomenon happen in disks or cores around high- and low-mass protostar(s) due to the large 
optical depths and higher temperatures only found close to the young star. Up to date detailed models of the 
ensuing chemical evolution at these much higher temperatures do not at present exist. 

It still is not clear why we do see an inner cavity on the disk in some molecules with transitions 
characterized by low and moderate excitation temperatures. This cavity maybe is not a real hole on the disk 
rather this could be due to photo-dissociation or opacity effects of these molecules toward the innermost parts 
of the disk. However, we do not discard the possibility that the cavity could be real and formed possibly by 
the tidal effects of the young multiple massive protostars in the middle of the "ring" as in the case of Ori 
139-409 located in Orion South where the protostars have lower masses (jzapata et al. 201ol ). The molecular 
emission with transitions characterized by high excitation temperatures in the lower energy states may arise 
from the possible compact circumstellar disks inside of the cavity that with the actual angular resolution 
are not resolved. The presence of a binary system at the center of W5 1 North has bee n suggested by the 
precession of the molecular outflow at very small scales found by Eisner et al.1 ((2002). This precession 
phenomenon may also explain the difference between the position angles of the blueshifted and redshifted 
components of the CO(7 = 3-2) outflow (Figure 1). 

The accreting disk hypothesis is supported by the detection of the powerful outflow that emanates 
with an almost perpendicular orientation to the object. Th e large disk s hows to be Keplerian and with a 
contracting velocity of about a few km s -1 as revealed by Eapata et al.1 (|2009T ). The models presented in 
Figure 4 were obtained using a slightly d ifferent P. A. = 30°, inclination angle of i=25°, a smaller size (3000 
AU), and without an inner cavity than in lZapata et al.1 (120091) . 



In Zapata et al. 



(2009) and here, we have modeled a contracting flattened disk in Keplerian rotation with 
a central hole, using the disk parametrization from Guilloteau, Dutrey & Simon (1999). The contraction is 
assumed to have the functional form of free-fall (i.e. V,„j oc -^=), with a reference velocity at the reference 
radius. However, our resolution is insufficient to determine the exact functional form. The model is for 
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the same molecules and transitions. A better fit than all other trials in our recurrence was found until we 
obtained similar structures in our model to those imaged (Figure 4). Most of the physical parameters in the 
model were constrained in the process. The model fits the observations reasonably well. 

In particular, the molecules HCOCH20H(62i349-62i2,5o) and HC3N(38-37) (v f =0) revealed the pos- 
sible presence of a warm companion located to the northeast of the disk (Figure 2). We suggest that the 
companion could be a consequence of disk fragmentation due gravitational instab i lities because the disk i s 
extremely large, see for a reference of this phenomenon iKratter & Matznen (l2006h : lKrumholz et al.l (12009). 



4.2. A possible sequence on the evolution of nascent massive stars? 

The multilayer disk imaged in this paper reveals a possible link between the Hyper/Ultra-compact HII 
regions and the hot cores. These observations suggest a possible sequence on the formation of massive stars. 
We present this simple sequence in Figure 7, and discuss it as follows: 



Phase I: In the first phase, a large and massive pseudo-disk is formed together with a bipolar outflow 
from a large and dense core. The pseudo-disk is surrounded by an infalling envelope with accretion 
rates on the order of 10~ 3 M Q yr _1 . The molecular emission arising from the circumstellar pseudo- 
disk could be observed as a single structure without a well defined temperature gradient across it. The 
pseudo-disks observed toward these objects are not classical cicumstellar disks associated with low- 
mass stars which are centrifugally supported, rather these might be contracting pseudo-disks (see for 
example the large and massive contracting disks associated with t he high-mass protostars AFGL490 
and NGC 7538S: iNakamura et aDll99ll : 1sandell & Wright! boioh . The pseudo-disk could be also 
likely circumbinary or circum-multiple because of the large multiplicity presented on the massive 
stars. 

Phase II: The large and massive pseudo-disk formed earlier is observed now in molecular shells 
or layers. This object is still surrounded by an infalling envelope. A clear temperature gradient is 
observed across the disk. The inner cavity could be formed by photodissociation due to the high 
temperatures close to massive stars or maybe to opacity effects. In this phase, no free-free emission 
from a HCHII region is detected probably because the infalling molecular gas partially quenches its 
formation or maybe because the massive central protostar(s) are/is in an early phase with not high 
enough temperatures to ionize their/its surroundings. W51 North is found here. Another possibility 
is that maybe there is cluster of late B-type in middle able of forming a strong multilayer molecular 
disk, but not for developing HCHII regions. 

Phase III: A HCHII region appears in the middle of the large multilayer pseudo-disk. At this stage 
the central massive protostar has the sufficient hig h temperatures (> 10 4 K) to ionize its surroundings, 
see fo r example the cases of NGC7538 IRS 1 (Franco-Hernandez & Rodriguez! |2004|; ISandell et al. 
20090 and MWC349 (ITafoya et al.ll20040 . These sources are maybe variable on time because of the 
strong accretion from the large disk continues toward the massive central star, see also G24.78+0.08 
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Al (balvan-Madrid et alJliooi) . There are some sources that might be situated at this phase, see for 
example, th e molecular rotating torpids obser ved around the UCHII reg i ons G 24.78+0.08 Al and 
G5.89-0.39 dSu et al.ll2Q09l : beltran et alJbood) . Moreover. Elaassen et all J2009h found in a group of 
HCHII and UCHII regions rotating molecular structures surrounding them. Therefore, the molecular 
toroids or rings around an UCHII are likely the same structures as in Phase II, but just hotter. The 
holes observed in these toroids might be due to the photodissociation of some molecules, or maybe to 
the formation of cavities by the presence of a multiple system of massive stars as mentioned before. In 
thi s phase the masive yo ung star(s) might be still accreting ionized material as proposed and described 
bv lKeto& WoodlJ2006h . 



Phase IV: In this phase the HCHII region expands and ionizes the remnant molecular disk or toroid, 
and creates an extended H II region surrounding one or maybe m ultiple young massive sta r s. Th e 
sequence from UCH II to extended HII regions has been discussed by lGarcia-Segura & Franco ( 2004 ): 



Franco et al. 



(2007). If the central star(s) is/are still accreting material, probably ionized material, an 



outflow show also be present even at this stage. 



5. Summary 

We have observed the massive and young protostar W51 North with the SMA using a high angular 
resolution and sensitivity at millimeter and submillimeter wavelengths. We give a summary of the results as 
follows: 



We report the presence of a large and single molecular disk around the object W5 1 north, and confirm 
the existence of a single powerful outflow emanating from the disk; 

The molecular emission from the large disk is observed in layers with transitions characterized by 
high excitation temperatures in the lower states (up to 15 12 K) being concentrated closer to the central 
massive protostar. The molecular emission at low or moderate excitation temperatures additionally 
exhibits a central cavity with an angular size around 0.7". This multil ayer disk connects the molec- 



ular ring and dusty compact disk reported recently toward this object (|Zapata et all 120091) . The LTE 



modeling presented in this study also confirms that the hot gas is only found in the innermost disk. 

A detail study of the kinematics of the molecular gas together with a model of a circumstellar disk in 
Local Thermodynamic Equilibrium shows that the disk is Keplerian, and is accreting fresh gas to the 
protostar possibly quenching the formation of a HCHII region; 

The thermal emission of the HCOCH2OH reveals the possible presence of a warm companion located 
to the northeast of the disk, perhaps produced by disk fragmentation or maybe is the interaction of 
the outflow a high density zone. The molecular emission of the SO and SO2 is observed in the 
circumstellar disk as well as in the outflow; 
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• We modeled all lines with a LTE synthetic spectra. For an assumed source size of 0.5 - 1 arcsec, 
our modeling yields a column density of 2 - 6 x 10 16 ~ 18 cm" 2 , a temperature of 250 - 800 K, and a 
linewidth of ~ 5 - 10 km s _1 . 

We thank the anonymous referee for many valuable suggestions. This research has made extensive use 
of the SIMBAD database, operated at CDS, Strasbourg, France, and NASA's Astrophysics Data System. 
This research made use of the myXCLASS program (https://www.astro.uni-koeln.de/projects/schilke/XCLASS), 
which accesses the CDMS (http://www.cdms.de) and JPL (http://spec.jpl.nasa.gov ) molecular data bases. 
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Fig. 1. Overlay of the SO2(222,20-22i_2i ) moment zero emission (color image), the 7 mm continuum emission from the compact dusty disk 

(brown contours ), and the CO(J = 3 - 2) moment zero emission from the bipolar outflow (blue and red contours). The SO2 and the 7 mm images 
were taken from lZapata et alj (2009). The integrated velocity range for the CO(J = 3 - 2) is from +20 to +50 km s for the blueshifted gas, and 
from +70 to +120 km s _1 for the redshifted gas. The brown contours are the 10%, 13%, and 16% of the peak of the 7 mm continuum emission, 
the peak is 100 mJy beam -1 . The synthesized beam of the CO(J = 3-2) contour image is 0.71" X 0.57" with a P.A. of 87.6° and is shown in the 
bottom right corner. The green triangle marks the position of the possible companion reported in this paper. The white dashed line shows how the 
blueshifted side of the outflow seems to be more likely energized from the central massive protostar W5 1 North rather than the possible companion. 
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Fig. 2. Spectra from the selected millimeter and submillimeter lines observed toward W51 North. The red line represents in all panels the 

synthetic spectrum obtained for the best fitting solution with the parameters shown in each panel. The spectra was obtained from the average of the 
total emission in each line using the task "imspect" of MIRIAD. 
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Fig. 3. Composite images of the molecular emission from the circumstellar large disk around the massive protostar W51 North. Left: The blue 

color scale and the contours are representing the integrated emission of the SO(89-7s), green ones the HC3N(38-37) (v f =0), and the red ones the 
HCOCH20H(62i3 49-62i2_5o)- The blue contours are from 5% to 90% with steps of 10% of the peak of the line emission, the peak is 53 Jy beam -1 
km s _1 , The green contours are from 30% to 90% with steps of 10% of the peak of the line emission, the peak is 26 Jy beairr' km s _1 . The red 
contours are from 30% to 90% with steps of 10% of the peak of the line emission, the peak is 3.5 Jy beam -1 km s . The synthesized beams of the 
millimeter and submillimeter observations are shown on the right bottom corner. Right: The light blue color scale and the contours are representing 
the integrated emission of the S02(19iig-18o ig), green ones the C2H50H(37g, 29-369,28), and the orange ones the HCOCH2OH(68i949-6720,4s)- 
The light blue contours are from 30% to 90% with steps of 10% of the peak of the line emission, the peak is 38 Jy beam -1 km s . The green 
contours are from 30% to 90% with steps of 10% of the peak of the line emission, the peak is 2 Jy beairr 1 km s . The orange contours are from 
30% to 90% with steps of 10% of the peak of the line emission, the peak is 27 Jy beam -1 km s -1 . The synthesized beam of the three measurements is 
the same and is shown on the right bottom corner. The triangles represent the position of the two possible embedded massive protostars as suggested 
by the hot molecular emission of the HCOCH20H(62[3 .49-6212,50) and HC3N(38-37) (v,=0). The white dashed line in both panels represents the 
position and orientation where the PV-diagrams shown in Figures 4 and 5 were computed. 
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Fig. 4. Position-velocity diagrams of the HC3N(38-37) (green scale color), and the HCOCH20H(68 19 49-6720,48) (pink scale color) line 

emission from the circumstellar disk. The PV diagrams were computed at a P.A. = 30° (see Figure 3). These diagrams are additionally overlaid 
with the position-velocity diagram at same P.A. of our Keplerian disk model (contours). The contours are from 20% to 90% with steps of 10% of 
the peak of the line emission of our model. The units of the horizontal axis are in arcseconds. The systemic LSR radial velocity of the ambient 
molecular cloud is about 60 km s -1 , see Table 1. The angular and spectral resolutions are shown in each panel on the bottom left corner. The 
synthesized beam of both images is 0.71" X 0.57" with a P.A. of 87.6°. The spectral resolution was smoothed to 1 km s - '. The origin in the 
horizontal axis is at R.A. = 19h23m40.05s, decl.= 14°31'05.0" (J2000.0). 
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Fig. 5. Position-velocity diagrams for the SO(89-7s) and the S02(19ij9-18o.i8) from the large circumstellar disk and the outflow. They were 

computed at a P.A. of 160° (see Figure 3). The contours in both panels are from 10% to 90% with steps of 10% of the peak of the line emission. 
The units of the horizontal axis are in arcseconds. The systemic LSR radial velocity of the ambient molecular cloud is about 60 km s~ 1 , see Table 1 . 
The angular and spectral resolutions are shown in each panel on the bottom left corner. The synthesized beam of both images is 0.71" X 0.57" with 
a P.A. of 87.6°. The spectral resolution was smoothed to 1 km s _1 . The origin in the horizontal axis is at R.A. = 19h23m40.05s, decl.=14°31'05.0" 
(J2000.0). 
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Fig. 6. SEP for the source W51 North combining 3.6 cm, 1.3 cm, 7 mm, 1.3 mm and 0.8 mm continuum data. The millimeter data were 

obtained from IZapata et al J |2009|) and the submillimeter data was obtained from Tang et al. (in prep.). The respective error bars were smaller 
than the squares an d are not presented here. The line is a least-squares power-law fit (of the form S v oc v a ) to the spectrum. Note that the value of 
a=2.8 obtained bv lZapata et alj|2009l) fit very well all millimeter and submillimeter measurements. The 1.3 cm point is a upper limit (5< x=5 mjy) 
obtained from Figure 1 of Eisner et al. 1 2002). The 3.6 cm point is a upper limit (5<x=4 mjy) obtained from Figure 5 of Mehri ngeJ jl994h . 
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Fig. 7. A very simple diagram showing the different evolutionary phases of a massive protostar. In this diagram, we do not included 

multiplicity for simplicity. 



